Adiabaticity and color mixing in tetraquark spectroscopy 
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We revisit the role ol color mixing in the quark-model calculation of tetraquark states, and 
compare simple pairwise potentials to more elaborate string models with three- and four-body 
forces. We attempt to disentangle the improved dynamics of confinement from the approximations 
made in the treatment of the internal color degrees of freedom. 
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I. INTRODUCTION 

There is a persisting interest in the quark dynam- 
ics applied to multiquark spectroscopy. The question 
is whether there exist compact hadron states beyond 
ordinary mesons (quark-antiquark) and baryons (three 
quarks) . 

In simple constituent models, several mechanisms have 
been proposed for binding multiquarks: chiral dynamics 
for light quarks, chromomagnetism, clustering of heavy 
quarks in a chromoelectric potential, etc. In this arti- 
cle, we concentrate on this latter effect, i.e., multiquark 
binding in a spin- and flavor-independent potential, and 
discuss the role of the internal color degrees of freedom. 
The validity of the existing models is beyond the scope 
of this note. In particular, we shall not discuss the tran- 
sition from color as a local gauge invariance to color as a 
global property of the wave-function in constituent mod- 
els. Still, even in its simplified version, color is a delicate 
ingredient of multiquark dynamics. 

The first studies on multiquarks within constituent 
models were based on simple color-additive potentials, 
extrapolated from meson and baryon spectroscopy. Al- 
ready for baryons, one hardly justifies the choice of a pair- 
wise interaction, half as strong as the quark-antiquark 
potential, but a more elaborate modeling, based on a 
connected F-shape flux tube linking the three quarks, 
does not change the results significantly. 

We shall follow in this paper this picture of a min- 
imal string linking the quarks. There are alternative 
non-trivial pictures of confinement, in particular the ones 
based on diquarks, which have been extended from the 
baryon sector to the multiquarks. See, for instance, [TH3]- 

The F-shape potential of baryons has been extended 
to tetraquarks and higher multiquark configurations [4]. 
There are multi-F connected diagrams in which the 
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string interaction links all quarks and antiquarks as a 
Steiner tree whose cumulated length is minimized. But 
the dynamics is dominated by the so-called "flip-flop" 
diagrams, with disconnected flux tubes for each of the 
quark-antiquark, three-quark or three-antiquark sub- 
clusters: the attraction comes from the minimum taken 
over all possible permutations of the quarks and of the 
antiquarks. 

This flip-flop interaction contains 3-body, 4-body, and 
higher-order terms, and is thus more delicate to handle 
in variational calculations. Moreover, when two quarks 
are exchanged, the color wave function is modified. In 
the latest studies 5-7 , this effect is not treated rig- 
orously. Instead, a type of adiabatic approximation is 
used: for any set of coordinates for the quarks, the po- 
tential is taken as the minimum of all permutations of 
the quarks and antiquarks, irrespective of the color wave 
function, and this minimum, as a function of the coordi- 
nates, is interpreted as an effective potential leading to 
a few-body spectral problem in which color has disap- 
peared. Interestingly, this strategy leads to stable multi- 
quarks for a large variety of constituent masses. This is at 
variance with the color-additive model, which binds only 
tetraquarks for large values of the quark-to-antiquark 
mass ratio. 

The question is thus whether the multiquark binding 
obtained in string models survives a non-adiabatic treat- 
ment of the internal color degrees of freedom. The prob- 
lem is analyzed in the present paper by reformulating the 
string-based interaction as an operator in color space. 

The outline is the following. In Sec. [TTJ we set the nota- 
tion for the color components of tetraquarks. In Sec. |III[ 
we present various models of tetraquark confinement with 
coupled channels in color space or an adiabatic approx- 
imation in which color disappears from the final wave- 
equation. The results are presented in Sec. |IV[ and the 
conclusions in SecfVl 
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II. COLOR STATES 

The internal color structure of tetraquark states is de- 
scribed in several papers, see, e.g. [8HIU]. We shall bor- 
row the notation of [5], in particular the names "true 
baryonium" (T) and "mock baryonium" (M), though the 
physics context is rather different in the present heavy- 
quark spectroscopy as compared to the color chemistry 
of the late 70s. 

The wave function for (1, 2, 3, 4) = (qqqq) is written as 

* = i' T \T) + IpM |M) (1) 

where |T) denotes a color state with the two quarks in 
a color 3 state, and the antiquarks in a color 3, while 
| M) corresponds to a color sextet in the quark sector 
and antisextet in the antiquark one. There is also the 
possibility of building the global color state out of quark- 
antiquark clusters either in a color singlet or octet. More 
precisely, we introduce 

\T) = |(12) g (34) 3 > , |M) = |(12)6 (34) 8 > , 

|1) = |(13) 1 (24) 1 ), |8) = |(13) 8 (24) 8 ), (2) 

110 = 1(14)1(23)!) , |80 = |(14) 8 (23) 8 ) . 

The relations between the different sets can be deduced 
from 




|iO = -^|T> + yfiM ); 

18')^ VI |T) + ^ |M) - 

Accordingly, the matrix elements of the potential in 
any basis are related to the ones in another basis, for 
instance, 

V u = (1|V|1) = ^Vtt+* V MM + Vtm , (4) 

and many similar relations. 

As stressed by Lipkin [TT] , in the limit of a tetraquark 
with two units of heavy flavor, (QQqq) with a large 
quark-to antiquark mass ratio M/m, the ground state is 
an almost pure \T) state, with the two flavored quarks in 
an antitriplet state, as in ordinary (QQq) baryons, and 
the two antiquarks neutralizing that color, as in (Qqq) 
antibaryons. In other words, the tetraquark state in the 
large M/m limit just uses well probed color structures, 
such as the 3 <g> 3 — > 3 coupling of two quarks in baryons. 

On the other hand, for smaller values of the mass ra- 
tio M/m, the simple models give at best a very shallow 
binding. Then the mixing of |T) and \M) is crucial to 
establish the stability. See, e.g., Brink and Stancu |12j . 



III. MODELS OF TETRAQUARK 
CONFINEMENT 

In early days of multiquark calculations, the potential 
was assumed to be pairwise, with the color dependence 
associated with the exchange of a color octet, namely 

V = -JeY,*i-hv(r ij ). (5) 

i<j 

Here, is the color operator for the i th quark, and is 
suitably modified for an antiquark belonging to the 3 
representation of SU(3). The normalization is such that 
v(rij) is the central part of the quarkonium potential. 

This model, however crude, gives at least the possibil- 
ity of studying the role of the internal color degrees of 
freedom inside a multiquark state. 

The additive model being subject to heavy criti- 
cism, an alternative was sought, inspired by the strong- 
coupling regime of QCD. It is referred to as the flip-flop 
model. For mesons, the potential is V — or, where r is 
the quark-antiquark distance, and a the string tension, 
which can be set to a = 1 by rescaling, without loss of 
generality. For baryons, this is the Y-shape interaction 
Vy = min a (ri a +7'2 a + r3 a ). For a tetraquark, the poten- 
tial is the minimum of the flip-flop term and a connected 
string, namely [H [5] 

V 4 = min(V s , V s ) , 

Vg = min(ri 3 + r 24 , ru + r 2 a) , (6) 
V a = min(ri a + r 2a + r ab + r 3b + r ib ) , 

a,b 

as pictured in Figs. [T| 
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FIG. 1. Top: schematic picture of the quark-antiquark and 
three-quark confinement. Bottom: three contributions to the 
tetraquark potential; the simple string model takes the mini- 
mum of these three contributions 

For each set of quark coordinates r iy the minimiza- 
tion in ^ implies a rotation in color space. This means 
that V4 is an adiabatic approximation, which tends to 
overbind the system. More importantly, the minimiza- 
tion, at least as it was carried out in [5] , does not account 
for any antisymmetrization. It just holds for distinguish- 
able quarks or antiquarks. 

As the model Q gives an interesting spectrum of sta- 
ble tetraquarks [5 , and, if extended to higher configura- 
tions, a spectrum of bound pentaquarks and hexaquarks 
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[51 [7] , it is crucial to estimate the amount of overbind- 
ing due to the adiabatic approximation, and the changes 
occurring when a proper antisymmctrization is imple- 
mented. 

We aim at constructing an operator in color space that 
tends to ^ in the adiabatic limit, at least when one of 
the three terms is clearly the minimum. For instance, 
if the (1,3) pair is clustered and lies far from the (2,4) 
pair which is also clustered, the potential is more easily 
described in the {|1), |8)} basis. However, working solely 
in this latter basis would require the choice of a value for 
the string tension between color-octet objects (there are 
studies within QCD, see, e.g., [13]) and an ansatz for the 
transition potential Vig. Instead, we shall combine the 
pieces of information coming from the singlet-singlet and 
triplet-antitriplet states, to deduce the full 2x2 matrix 
of the potential in the {\T}, \M}} basis, in which the 
four-body problem will be solved. 

More precisely, we will consider four different models: 

• Model A is the adiabatic limit given by ([6]), al- 
ready used in [5]. However, the Steiner-tree with 
two junctions, that plays a marginal role, is ne- 
glected. Hence, this is the pure flip-flop model. 

• Model B is a smooth version of the adiabatic ap- 
proximation. We use g(x) = 1/(1 + x n ) and its 
complement g'(x) = 1 — g(x), with some large but 
finite exponent n = 5 to soften the transition be- 
tween different limiting regimes of the string model. 
In practice, we replace the above Vtt by 



Vtt = 03 3 v tt + (1 - 03 9a) 



3Vvi> + Vu SVn+Vyv , 
: 9i + - A Si 



(7) 



where 

Vii = 7*13 + r 2 4 , Vvv = ru + r 23 , (8) 

and 



^33 \ , ( ^33 



0i = 1 - 9i = 9 



Vii 
V V1 . 



(9) 



Then, the potential is known in any basis from Vn, 
Vyy and Vtt- For instance, in the {|1), |8)} basis, 
one uses V\\ and 



^16 



1 



4^2 



{Vii + Wvv ~ 4V TT } 



(10) 



Vss = ^{-V ll + 3Vyy+2V TT } , 



and it is readily checked that if gi — > 1, the poten- 
tial becomes diagonal. The four-body problem is 



more conveniently solved in the {\T), \M)} basis, 
and besides Vtt, the relevant matrix elements are 



V T M = ^{Vn-Vyy} 



(11) 



Vmm — ^{3^ii 



• Model C is the color- additive model of normal- 
ized to a unit string tension for quarkonium. 

• Model D is the crude adiabatic limit of the previous 
one. This means that for any given set of positions, 
the 2 x 2 matrix consisting of Vtt , Vtm and Vmm of 
model 4 is diagonalized, and the lowest eigenvalue 
is taken as the effective four-body potential, irre- 
spective of any symmetrization or antisymmetriza- 
tion. Of course, model D is more attractive than 
model C. 

• Several other variants have been envisaged, but 
abandoned as leading to collapses or inconsisten- 
cies. For instance, it is tempting to use relations 
similar to ( 11 ) to express Vtm and Vmm from Vii, 
Vi'i' and Vtt, as given by the string model. But 
then the Hamiltonian is not bounded below. 



IV. RESULTS 

The results are shown in Table Q] The aim is not to 
provide a benchmark of four-body calculations. What 
really matters, is how the ground-state energy evolves 
when going from a model to another one. The variational 
estimate has been carried out using just a few Gaussians, 
and in the case of models B and C, imposing the proper 
symmetry. 

For this purpose, we make use of a wave function with 
the relevant symmetry for each color component. As the 
color vector \T) (\M)) is antisymmetric (symmetric) un- 
der both the exchange of the identical quarks and the 
identical antiquarks, it has to be combined with a ra- 
dial wave function with the proper symmetries. The way 
of constructing these wave functions has been explicitly 
detailed in Ref. [TO]. We will just draft here its main 
characteristics. The radial wave function is taken as a 
linear combination of generalized Gaussians depending 
on six variational parameters a%j of the form, 



4 

E 

k=l 



a k exp 



3 

i> 3 =i 



(12) 

where s k are six-component (ij) vectors made of arrange- 
ments of positive and negative signs. Once combined 
with the proper election of the signs ak, they give rise 
to radial wave functions with the following symmetries 
in the radial space under the exchange of quarks and 
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TABLE I. Results for the various models A, B, C and D, as 
a function of the quark-to-antiquark mass ratio M/m 



M/m 


A 


B 


C 


D 


Threshold 


1. 


4.644 


4.803 


4.702 


4.596 


4.676 


2. 


4.211 


4.306 


4.275 


4.160 


4.248 


3. 


4.037 


4.131 


4.112 


3.984 


4.086 


4. 


3.941 


4.041 


4.010 


3.891 


3.998 


5. 


3.880 


3.985 


3.954 


3.828 


3.942 


10. 


3.742 


3.860 


3.834 


3.685 


3.831 



antiquarks: SS, SA, AS, or A A, where S stands for sym- 
metric and A for antisymmetric. 

The results from models A and D are very similar for 
the variational energies. This means that the extra at- 
traction noticed in [5] is mainly due to the color mixing in 
the adiabatic approximation. Restoring the interaction 
as an operator in color space is less favorable for multi- 
quark binding, and it is readily seen that models B and 
C led to comparable predictions. The color structure of 
the wave-functions is also rather similar in models B and 
C, with mostly a singlet-singlet configuration. 

However, if one looks at the details, it can be realized 
that the flip-flop and and the additive models differ. In 
particular, the average QQ separation, (fia) is smaller 
in the additive model than in the flip-flop model. If the 
mass ratio M/m — > oo, this effect becomes more pro- 
nounced, with (1*12) — > in the former case and (7*12) 
evolving to a finite value (for m fixed) in the latter one. 
This influences the contribution of the VF-exchange con- 
tribution to the weak decay of (bcqq), which also plays a 
role in the decay of (bcu) baryons |14) . 

This large M/m limit is rather interesting. A detailed 
numerical study would involve the computation of an ef- 
fective QQ interaction, V e e(ri2), sum of the direct QQ 
interaction and of the qq binding energy around them, 
very similar to the Heitler-London potential of two pro- 
tons in the adiabatic treatment of the hydrogen molecule. 
Already, for baryons with two heavy quarks, (QQq), it 
was stressed that the Born-Oppenheimer approximation 
works very well to reproduce the properties of the lowest 
states [IS]. For (QQqq), it is striking that the effective 
potentials I4ff(ri 2 ) have qualitative differences in the ad- 
ditive model and in the flip-flop one. In the additive 
model, the minimum of T4ff(ri2) is reached at r*i2 = 0. 
The qq contribution is stationary there. So, up to a con- 
stant, the potential is dominated by the direct QQ term, 
which is r 12 /2 here. Thus the average separation de- 
creases as M^ 1 / 3 , according to the well-known scaling 
laws in a linear potential |16j . In the flip-flop model, 
Ves(fi2) is minimum at some finite distance. This can 
be seen directly. This is also compulsory, if one wishes 
to understand why (QQqq) is stable in the large M/m 
limit, as seen by the following reductio ad absurdum: 
suppose that in the Born-Oppenheimer limit of the flip- 
flop model, (r 12 ) — > 0, then the flip-flop energy of (QQqq) 
and the energy of its threshold, (Qq) + (Qq), which are 




FIG. 2. Flip-flop potential for ri2 — > and cumulated poten- 
tial of two mesons 

pictured in Fig. [2] and given by 

V c ff(r 12 ) = min(n 3 + r 24 , ru + r 23 ) , 
Vth = ri3 + r 24 , 

would coincide if r\i — and the inequality among ener- 
gies (QQqq) < (Qq) + (Qq) would be impossible. 

V. CONCLUSIONS 

Our results illustrate the role of antisymmetrization 
in preventing from a proliferation of multiquarks. In 
particular, the main difference between the naive color- 
additive model and the crude flip-flop model comes from 
the treatment of color. It remains that in such a flavor- 
independent confining, the binding of (QQqq) is obtained 
for heavy enough quarks. 

More ambitious tools are obviously needed to handle 
the four-quark dynamics. In particular, in the double- 
charm sector, the binding effect of the confining interac- 
tion is probably not sufficient. In addition, the spin-part 
of the wave function can be made favorable: as stressed, 
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FIG. 3. Some weak decays of the (bcu) baryon and (be) meson 
leading to a tetraquark in the final state 
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e.g., in [T7] and references therein, there is a light -light 
interaction in (ccud) which is absent in the two mesons 
constituting the threshold. 

The production of such tetraquarks can be accessible in 
B-factories [T7] and at the proton-proton colliders. Note 
also that (ccqq) can be a decay product of hadrons con- 
taining charm and beauty. From a (bcu) baryon, for in- 
stance, the Cabibbo allowed b —> c + W~ —> c + d + u, 
combined to a dd pair creation, leads to (bcu) — > (ccud) + 
(ddu). See Fig. [3] From (6c), one could first envisage the 
Cabbibo suppressed b — > u + W~ — > u + c + s, and after 
the creation of a light quark-antiquark pair, this would 
monitor a decay (be) — > (ccud) + (sd). Of course, the 
CKM suppression factor is rather effective here. Perhaps 
more promising is the chain b — ¥ c + W~ — > c + c + d 
giving altogether (ccuudc) after a uu pair creation. This 
could lead to D~ +X, where D is an anticharmed meson 
and X one of the new hidden charm resonances reviewed, 
e.g., in [T5], [TH] ■ Another combination is (ccud) + (cu), 



with, however, a different topology of the quark diagram 
and thus different color and OZI suppression factors, as 
discussed by Lipkin in a different context [20J . Anyhow, 
any heavy-quark factory should lead to the discovery of 
heavy tetraquarks with suitable triggers. 
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